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. . . We calculate the self-energy and the vertex radiative corrections to the effect of parity noncon- 

\ servation in heavy atoms. The sum of the corrections is of the form ^ln(Ac/''o) + B, where A 

. and B are functions of Za, and Ac and ro are the Compton wavelength and the nuclear radius, 

' respectively. The function A is calculated exactly in Za and the function B is calculated in the 

leading order. In the leading order A oc a(Za)^ and B oc a[Za). The sum of the corrections is 
^ , —0.85% for Cs and —1.48% for Tl. Using these results we have performed analysis of the experi- 

^ • mental data on atomic parity nonconservation. The values obtained for the nuclear weak charge, 

Q \ Qw = —72.81{28)exp{36)theor for Cs and Qw = — 116.8(1.2)e2,p(3.4)theor for Tl, agree with predic- 

tions of the standard model within 0.6a. As an application of our approach we have also calculated 
\l ■ dependence of the Lamb shift on the finite nuclear size. 
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I. INTRODUCTION 



■ Atomic parity nonconservation (PNC) has now been measured in bismuth [y , lead , thallium , and cesium . 
] Analysis of the data provides an important test of the standard electroweak model and imposes constraints on new 
'"q^ physics beyond the model, see Ref. The analysis is based on the atomic many-body calculations for Tl, Pb, and 
I _ Bi § and for Cs (?[ |] (see also more recent Refs. 0)- Both the experimental and the theoretical accuracy is 
Oh. best for Cs. Therefore, this atom provides the most important information on the standard model in the low-energy 
sector. 

In the many-body calculations ^ the Coulomb interaction between electrons was taken into account, while 

the magnetic interaction was neglected. The contribution of the magnetic (Breit) electron-electron interaction was 



calculated in the papers [|ll|, 12 1. It proved to be much larger than a naive estimate, and it shifted the theoretical 
prediction for PNC in Cs. 

■ Radiative correction to the nuclear weak charge due to renormalization from the scale of the W-boson mass down 
to zero momentum has been calculated long time ago, see Refs. |lj]. This correction is always included in the 
analysis of data. However, another important class of radiative corrections was omitted in the analysis of atomic 
PNC. This fact has been pointed out in work that demonstrated that there are corrections ^ Za^ caused by the 
collective electric field of the nucleus. Here Z is the nuclear charge and a is the fine structure constant. The simplest 
correction of this type is due to the Uehling potential. It has been calculated numerically in Ref. ||l^ and analytically 
in our paper p^ . In that paper p7| we have also analyzed the general structure of the radiative corrections caused 
by the collective electric field. It has been shown that, as well as the usual perturbative parameter Za, there is an 
additional parameter ln(Ac/f'o) where Ac is the electron Compton wavelength and tq is the nuclear radius. 

In the present work we consider radiative corrections to the atomic PNC effect due to electron self-energy and 
vertex. The total correction is of the form A\n{Xc/ro) + B, where A and B are functions of Za. We calculate the 
function A exactly in Za and the function B in the leading order. In the leading order A cx a{Za)'^ and B oc a{Za). 
Results of the leading order calculations have been reported in Ref. flq] . Using results of our calculations we reanalyze 
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the experimental data for atomic PNC. Agreement with the standard model is excellent. 

As an application of our approach we have also calculated the dependence of the Lamb shift on the finite nuclear 
size. Agreement of our analytical formula with results of previous computations |l9| p(| is perfect. Structure of the 



IV 



present paper is following. In Sec. || we discuss the general structure of the PNC amplitude. In Sees. Ill and ^ ^ 

we calculate logarithmic parts of the self-energy and the vertex contributions, respectively. In Sec.^ we calculate the 
linear in Za correction to the PN C am plitude. In Sec. VI we calculation the dependence of the Lamb shift on the 
finite nuclear size. Finally in Sec. VII we analyze experimental data on atomic PNC and present our conclusions. 



II. GENERAL STRUCTURE OF THE PNC AMPLITUDE 



The strong relativistic enhancement is a special property of the atom ic PNC effect. The relativistic enhancement 
factor is proportional to i? ~ {\c I ZarQY'^^~^^ ^ where 7 = -^/l [ZaY' . The factor is R «3 for Cs and i? « 9 for Tl, 
Pb, and Bi For nuclear radius we use the formula 

ro ~ 1.1 J^l^ fm w 1.5 Z^/^ fm. (1) 

The relativistic enhancement factor R is divergent at ro ^ 0. This makes the relativistic behavior of PNC very much 
different from the behavior of the hyperfine structure. The logarithmic enhancement of radiative corrections mentioned 
above is closely related to the existence of the factor R. The Feynman diagram for the leading contribution to the 
PNC matrix element between pi/2 and Si/2 states as well as diagrams for one loop radiative corrections are shown in 
Fig.|l|(a) and Fig|^(b-f), respectively. Strictly speaking the diagrams in Figs.|l|(b-f) are not quite well defined because 
they describe the matrix element between states with different energies, so it is not clear what energy corresponds to 
the external legs. However the external legs describe states of the external atomic electron, say 6s, 6p, 7s,... states in 
Cs, that have very small binding energies of the order of ma^ (m is the electron mass). Therefore, the uncertainty in 
the definition of the diagrams appears only in the order a^{Za) which we do not consider in the present work. The 
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FIG. 1: (a) Leading contribution to the PNC matrix element, (b-f) One loop radiative corrections. The double line is the exact 
electron Green's function in the Coulomb field of the nucleus, the cross denotes the nucleus, the zigzag and the dashed lines 
denote Z-boson and photon, respectively. 



diagram Fig.y(b) corresponds to a modification of the electron wave function because of the vacuum polarization. 
This correction calculated analytically in Ref. jl^ reads 

S,^a (^^Za + [ln'{bXc/ro) + /]) , (2) 

where b = exp(l/(27) — C — 5/6), C ~ 0.577 is the Euler constant, and / ^ 1 is some smooth function of Za 
independent of rg. Hereafter we denote by S the relative value of the correction. So, Eq. (^) represents the ratio of 
diagrams Fig.|l](b) and Fig.0(a). 

The renormalization of the nuclear weak charge Qw from the scale of the W-boson mass down to g = was 
performed in Refs. jl3[ However, as has been pointed out in Ref. [|l^, atomic experiments correspond to 

q I/tq 30MeV, see also Ref. [|2) . The correction due to renormalization from q = to q = 1/ro is described by 
diagrams Fig.0(c) and Fig.0(d). It has the form 

Scd = -^(l-4sin2 0v^)ln(Ac/ro) ~ -0.1%, (3) 

where 9w is the Weinberg angle, sin^ 9w ~ 0.2230, see Ref. 
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Diagrams Fig.0(e) and Fig.|](f) correspond to the contributions of the electron self-energy operator and the vertex 
operator, respectively. Neither of these diagrams is invariant with respect to the gauge transformation of the elec- 
tromagnetic field. However, their sum is gauge invariant. It has been demonstrated in Ref. ||l^ that the correction 
6e + Sf is of the form 

d,f^S, + Sf=AHbXc/ro)+B , (4) 

where A and B are functions of Za, and the constant b is defined after Eq.(^. In the leading approximation in the 
parameter Za the functions are, A = {0.2 /t^) a{Za)^ and B = ai a{Za). In the work we have also obtained 
preliminary estimates for the coefficients Oi and 02. The estimates were based on an assumption of analogy between 
the polarization operator and the self-energy operator. This assumption and hence the preliminary estimates proved to 
be wrong. In the present work we calculate the coefficients ai and 02 exactly. Moreover, we also calculate the function 
A exactly in the parameter Za. The functions A and B are calculated by different methods. For calculation of A we 
use electron Green's functions in the Coulomb field at r <C Ac and the Feynman gauge. The non-logarithmic term B 
we find using the effective operator method and the Fried- Yennie gauge . The method is based on calculation of 
the low-energy scattering amplitude in the potential induced by Z-boson exchange. 

Let us consider the PNC matrix element between pi/2 and Si/2 states. The wave function of the external electron 
is of the form 

where £7 and J7 = — (cr • n)f2 are spherical spinors[p4|. Calculation of the weak interaction matrix element in the 
leading approximation gives |^ 

< Pi/2|(2V2)-iGi.Q^p(r)75|si/2 >o= A/q a {F,Gp - G,Fp)|,=,„. (6) 

Here Gp is the Fermi constant, and p{r) is the density of the nucleus normalized as / p{r)dr — 1. At small distances, 
r <C ZaXc, the electron mass can be neglected compared to the Coulomb potential energy, and solution of the Dirac 
equation for the radial wave functions reads 

F = Nr-^-^, G = N-^r''-^^ (7) 

K — 7 

where k = —1 for Si/2-state and k = I for p]^/2"State, and N is some constant dependent on the wave- function behavior 
at large distances (r ~ as) Q. Calculation of the constant is a very complex many-body problem, see Refs.|Q |j. 
Fortunately, the constant is cancelled out from the relative correction which we consider. Due to Eq. the matrix 
element is divergent at tq ^ , Mq cx Tq^ ^. This leads to the strong relativistic enhancement factor mentioned 
above. 

The correction to Mq related to the diagram Fig|^(e) can be obtained as follows. The electron self-energy operator 
S being substituted to the Dirac equation, m — > m+S, leads to the Lamb shift of energy levels and to the modification 
of electron wave functions, see, e.g., Ref. The wave function modification infiuences the matrix element of the 

weak interaction. We will call the Dirac equation with account of the self-energy operator the modified Dirac equation. 

It is convenient to search for solution of the modified Dirac equation in the form (||) with substitutions F 
F{1 + G ^ G(l -I- G(i)), where F^^) and G^^) are corrections due to the self energy. Using this form in Eq.(|^) 

together with we obtain the following relative correction to the PNC matrix element 

S = m(i)/m(o) = ^(1) + G«*l + [g« + 41)*] 1 . (8) 



Similar to Eqs.(||),(^, we represent the self-energy term in the Dirac equation 'f^ — Hip as 



^'^^ L-^(2)^pl r''-' . (9) 

K — 7 

In the next Section we will show that F^^'' and G*^^^ are independent of r at r <C ZaXc- The functions F'^^'> and G'-^-' 
satisfy the following equations 

^^-F(^)+G«=-^G(^) 
J + K OX Za 

1 ^^^'^ j7(i)+GW = , (10) 



K — 7 dx Za 
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where x 



ln(r/Ac). The solution of Eqs. ( |To| ) reads 

Zax „/T\-, 1 



27 



[G(2) 



1 



[(«: + 7)G(2) 



- 7)i^(2)] _^ „ ^ 



(11) 



where a is some constant. Note that in the leading logarithmic approximation (the terms proportional to a;) F^^^ 
G(i). Substituting (0) into Eq.(|), we find 



ZaL 



G(,2)* 



const. 



(12) 



Here L = ln(feAc/ro), and by const we denote Za-dependent terms that do not contain the large logarithm L. It will 
be shown below that F, 



:2) 



G?^ and G 



(2) 
P 



Fs"^' . Hence, the correction 5^ is of the form 



be = Ke 

7 



G(2) _ f(2) 



L + const . 



(13) 



f2l (2) 

We will demonstrate that Gs and Fs are odd functions of Za. Therefore, the factor before the logarithm L is an 
even function of Za and hence the leading term of Za-expansion of this function is proportional to Q;(Za)^. The const 
in Eq. ( |l3|) can be determined from the condition that the correction to the wave function is orthogonal to the solution 
of the non-modified Dirac equation. Therefore, the main contribution to the const comes from the distances r ~ Ac, 
and hence the leading term of Za-expansion of the const is proportional to a(Za). The correction Sf corresponding 
to the diagram in Fig.|l](f) has the same properties as Se- 



III. LOGARITHMIC CONTRIBUTION OF THE SELF-ENERGY OPERATOR 

At r ^ Ac one can neglect the electron mass in the propagator and write down ^'^ = S-f/; as 

00 

«'^(r2) = 2a J de J dri^^g{r2,ri\ie)-f'' D{r2,ri\ie)^jj{ri) , 



(14) 



where Q is the electron Green's function in the Coulomb field, g^,jD is the photon Green's function in the Feynman 
gauge. In ( p^ ) we have deformed the contour of integration over the energy e in such a way that it coincides with the 
imaginary axis. According to the integral representation derived in ||2^, the Green's function is of the form 



C/(r2, ri| ie) = > / ds exp [ 2zZaAs — fc(ri + r2) coth s ] T , 

47rrir2 f-' J 



1=1 



T - 7"[1 - (7n2)(7«i)] [A|/L(y) - iZ a coih s h.{v) 
'k{r2 - ri) 



B + 7°[1 + (7n2)(7«i)]A/2.(y)A 



2 sinh^ s 
2ky/rrr2 



■(7,ni + n2)B + coths(7, n2 - ni)A 



smh s ax 



hu{y) , 
d 



B = -{Pi{x) - Pi^^{x)) 
dx 



(15) 



Here fc = |e| , A = e/fc , ■n,i_2 = »"i,2/''i,2, x ^ rii ■ n2 , Pi{x) is the Legendre polynomial, is the modified Bessel 
function of the first kind, v = yJP — (Za)^. Using relation 



J {ni X n2)^{ni ■ n2)dni ^ , 
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that is valid for arbitrary function <J>(a;), and Eqs.(|^), (|T5|), and(p^ we get for s-wave 

Fs^\ iar ^ ^ f f -2 f 

^« / i=ix=±i{ J { 

fs = il + x)B [\^lUy)~iZacothsl2Ay)] + (1 - x)AXl2Ay) 



iZa 
1 + 7 



K'T-i - ri) 
2 sinh'^ s 



(1 — a;)Acoth s 



hu{y) 



(1 + x)(l - 2x)B yX^I'^M - iZacoihsl2.{y)\ " {I - x){l + 2x)AXl2u{y) 



iZa 



[l + x)B 



k{i^2 - ri) 
2 sinh^ s 



+ (1 - a;)Acoths 



hAy) 



Performing similar calculations for p-wave we find that Fp 
photon Green's function D reads: 



(2) 



gI^^ and G 



(2) 



n(2) 



(16) 



The explicit form of the 



D(r2,ri\ie) = - 



,-kR 



R=\r2~ri\ 



(17) 



Using Eqs.(|l6|), (|T^),(^3|) as well as the substitution ri — r2P , e — E/r2 we obtain the following expression for the 
self-energy contribution to the function A (see, Eq. 



AsE = j'^^ J ^^^^ Jdsexp[-E{p+l)coths- ER] (Ai + A2 

^—1 n n 



Ai 



x{l + x)B |^|/2^(j/)(ZQ;)sin0 - [Zaf cos coth s 72,. (y) + {Za){l - x^)Aiiin(j) huiy) 



At — cos ( 



E-i(l - p) 
{l + x)B— — + (1 - 2;)Acoths 
2 sinh s 



h.{y). 



(18) 



Here (f) — 2Zas , R = yp^ — 2paH-l, j/ = 2Ey/p/ sinh s. It is obvious from this formula that Ase is an even function 
of Za. 

Equation ( [T^ ) is not quite correct. The point is that the integrand in this equation has been derived for ri 7^ r2. 
However, there are also terms containing i5-function 6{r2 — Ti) and derivative of this (5-function. To account for 
these terms we represent the self-energy operator as a series in powers of the Coulomb field of the nucleus, S = 
So -I- Si -I- S2 + see FigJ|. The first two terms in the series, Sq and Si, require the ultraviolet regularization. It 



X XX 

FIG. 2: The electron self energy expanded in powers of the Coulomb field. The solid line is the free electron Green's function, 
the cross denotes the nucleus, and the dashed line denotes the photon. 

means that they contain terms with the 5-function and derivative of the (5-function. Both Sg and Si depend on the 
parameter of the ultraviolet regularization. However, the contribution of Sqi = Sq -|- Si to Ase is independent of the 
regularization parameter because of the Ward identity. The operators S„ (n > 2) do not require any regularization. 
Note that Ai in Eq.(p^) corresponds to the contribution of the operators S2n+i, while A2 corresponds to that of the 
operators S2n. Let us represent the self energy operator as S = Sqi -|- S(">^). Corresponding contributions to the 
function Ase are A^g^ and A''g^^\ hence 

AsE-A"'J+A^->'\ (19) 

The terms singular at ri = r2 contribute only to A'g^ . To overcome the problem with singularity we calculate A'g^ 
using momentum representation and the standard covariant regularization of Sq and Si . A simple calculation in 
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the momentum space gives the following expression for Sq 



So(P2, Pi) = -TT- 
ZTT 



1 



A 



(20) 



where A is the parameter of the regularization, and 'd{x) is the step-function. Because of the Dirac equation, pij} 
(—Za/r)"fQ^, the operator Sq in the modified Dirac equation can be replaced by 



So(P2, Pi) 



2a{Za) 

{P2~Plf 



1 



A' 



7 + TTln(^M(A^-P^) 



(21) 



The operator ( pl[ ) has the same 7-matrix structure as that of Si. As a result the sum of Si and So is independent 
of A : 



Soi(i>2, Pi) 



a{Za) 

(P2 -Pl) 



1 1 



1 + 2 



dxdy 



[x^y^{P2-Pif + il~x){xpl+pip2)] >7" , 



T = y(l - xy)pl + {1 - y){l - X + xy)pl ~ 2xy{l ~ y)pi ■ p2 



(22) 



The Fourier transform of Eqs. (^,(0) gives the following result for the Si/2 wave function in momentum representation 



Hp) = N 



47rr(7 + 1) 



cos(7r7/2)r2 
{Za/'^) sin(7r7/2)(cr • np)fl 



where rip = p/p- Using Eqs. 



Sl), and definition of F^^^ and G^^^ given in Eq. (ph, we find 



^(2) ^ cos(7r7/2) 
sin(7r7/2) 



P2 v. / N d-Pi 
^:p2^oi(P2,Pi) 

Pl 



(27r)3' 



(72 - l)sin(7r7/2) 
7COs(7r7/2) 



Pi 



pI / / N dpi 
:^(Pi • P2)Soi(P2,Pi) 



(27r)3 



Evaluation of these integrals and substitution of Fs and Gs in Eq. ( p^ ) (see also Eq 



(01) _ 
SE — 



a{Zaf 
27r7 



V'(2 + 7) + ^(2-7)-^( 



+ 7^(1/2)- V'(l)-l + ^ + 
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+ '^^ 

— )-^(l-2)- 



gives 

1-7 



2(1 + 7) 



7^ 1 — 7^ 7sin(7r7) 



(ilnr(a;)/da;. The first term of Za-expansion of the function A^g]^ is —a{Za)'^/'K. Note that even at 

(01) 



(23) 



(24) 



(25) 



90 the exact value of -4^^'' differs from that given by the first term of the expansion by 10% only. 



Here il}{x 
Z 

To find one has to subtract from (|18| ) the contribution (pS]). To do this we subtract from the integrand in 

( p^ ) the first term of its Za-expansion at fixed 7. The point is that 7 in the integrand comes from the wave function, 
but not from the self-energy operator. Therefore, doing this subtraction one has to consider 7 as an independent 
parameter. After the subtraction we make substitutions E — £/2y/p and p — exp(2T), and take the integral over the 
azimuth angle of the vector p. Altogether this gives 



A 



(n>l) 
SE 



dx 



de dr ds 
V 



exp [— e(I? + coshr coths)] A 



000 



A = 



x(l + a;)Bcosh(27r) {(Za)| [sin0/^,(y) - ^^(y)] - (Za)^ coths [cos0/2.(y) - hiiy)]} 

+ {Za){l - a;2)Acosh(27r) [sin0/2.(y) -^hi{y)] 

, £7 sinh T sinh(27r) 



[cos(j)l2v{y) - hiiy)] 



-{l + x)B- 



2 sinh s 



(1 — .■E)ylcoths cosh(27r) 



(26) 



7 



Here I? = y sinh^ r + (1 — a;)/2, and y — ej sinh s. The leading term of Za-cxpansion of J^g^^^ 

-Q(Za)2(7r2 -9)/(67r). 

It corresponds to the contribution of the operator S2- According to Eq. (|l^), the leading term of Za-expansion of 
AsE reads 



(-^ - 3) 



(27) 



Numerical calculation of integrals in Eq. (g^) together with Eq. ( |25D gives the function Ase exactly in Za. This 
fimction is plotted in Fig. 




FIG. 3: The junction Ase calculated in all orders in Za. Value of the function is given in units "^^^^ 



IV. LOGARITHMIC CONTRIBUTION OF THE VERTEX 

Let us consider now the logarithmic contribution of the vertex operator (diagram (f) in Fig.|l|). According to the 
rules of the diagram technique, the expression for the vertex correction to the matrix element is of the form 



Mf = 2a jde I I dridr2dr3'iljp{r2)^f,g{r2,r3\ie) 



X 7^7''p(r3)^ (ra , ri | ie)-/^'ips (r i )!) (rs , ri\ie) 



(28) 



where, as above, we neglect in (^8|) the electron mass. We use expansion of the electron Green's function in spherical 
waves, see Eq. (p^. Because of the contact nature of the weak interaction, only the term with angular momentum 
j — 1/2 is important in the expansion, this term corresponds to / = 1 in Eq. (|l5|). Note that j = 1/2 dominates 
only in the calculation with logarithmic accuracy. For a calculation of the constant near the logarithm one has to 
take into account all the partial waves in spite of the contact nature of the weak interaction. Now we calculate only 
the logarithmic contribution. The logarithmic contribution to the amplitude Mf comes from the region of integration 
ri ~ r2, Ac ^ ?'i,2 ^ ^ ro, € ^ ^/ti,2- In this case the parameter s in the integral representation Eq.(|l5|) is of 
the order of unity, s ~ 1 , and hence the argument of the Bessel function is small, y ~ ('"3/''i.2)^^^ ^ 1- Expanding 
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Bessel functions we transform Eq.(|T^ to the following form 

G{r2,ri\ie) = — — — — / ds exp[2iZaXs - k{r2 + n) coths] 

47rrir2r(27 + Ij Jq 

X (^7^) ^|[l-(7n2)(7«i)]7"(7A-i^acotlis) + [l + (7n2)(7ni)]7"A 

— —— — (7,n2 - (7,n2 -ni)cotlis L (29) 
2 smn s I 

It is convenient to perform integration by parts over s in the term proportional to (r2 — ri) . After that we substitute 
( p9| ) in equation ( |2^ ) and take the integral over ra. As a result we get 

„ , ^ , OO CXXJO 

ar f f f dridr2 > ff 

Mf = 2^2p2° 2^ 4, 1) y J J — X! yy c?si'^'S2 exp [2iZQ;A(si + S2) - e(ri cothsi + r2 coths2) ] 

A = ±l Q g 



^1^) ^p(r2)|7°[l - (7ni)(7«2)] (1 - cothsi coths2) +7°[1 + (7"i)(7W2)] 



\ sinh s 

X [7 (1 + cothsi coths2) — iZaX (cothsi + coths2)] + A(7, rii + 722) (cothsi — coths2 



+ (7, Til — 7x2) [7A (cothsi + coths2) — iZa (cothsi coths2 + 1)] j7^7/'s(ri) D{r2,ri\ie) . (30) 

Then we substitute the photon propagator (|l7| ) and the electron wave functions (||), (0), and find the following 
expression for the relative correction 5f 



u. 1,11 dridr2 expi—tR) ■r-^ / / m n 

= ' 87r3P(27+ 1) 7^^7 'Wl ^ ^ / ds,ds2e^p[2^ZaX{s, + S2) 



00 /'OO 







Jo 



/ eVir2 A^"" 

-e(ri cothsi + r2 coths2) ] 

\smhsi smhs2 / 

= (1 — a;^)(l — 7) [(1 + 7) cothsi coths2 ^(1^7) — iZaX (cothsi + coths2)] + 2(1 — xcothsi coths2) 



R = y rl — 2xrir2 , x = rii ■ n2 . (31) 

At the next step we introduce new variables p, t, E via the relations ri = pt, r2 = p/t^ ^\frW2 = -E-, perform 
integration over p with logarithmic accuracy in limits ro < p < Ac, and take integrals over all the angles except of 
the angle between vectors r\ and r2- This gives the following expression for the contribution of the vertex operator 
to the function A defined in Eq.(0): 



00 CXD 



Ay = - ^^2(27 + 1) /'^^ jl j ^""P*-^^^^ ^ // '^'^^'^'^^ exp[2^ZaA(sl + S2) 



00-1 -"^-^I'oo 

?2 \ 27 



-£;(icothsi + (l/<)coths2)] — ) J" ■ (32) 

\smhsi smhs2 / 

Here R — ^JW~^^TJW^^2x and J- is defined in (^l]). A simple integration over E leads to 

2ar(47+l) dt dx ^ f""/""" dsids2 exp[2iZaA(si + S2)] 
~ ~7rr2(27+l) io ~ 1-1 'R ^tiJo Jo (sinh si sinh 52)2-^ 



J" 



[t COth Si + {l/t) COth S2 + 



(33) 
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The last integration which can be performed analytically is the integration over x. The result is rather cumbersome: 

Qr(47 + 1) dt dsids2 e^pl'iiZaXlsi + S2)] 

T 2^ 



Av = 



7r7r2(27 + l)7o ^ 1^0 ^0 



(sinh si sinh 82)^'' 
(1 — cothsi coths2) (1 + cothsi coths2) 



1 



47-1 



\t - l/t\{Q + cothsi coths2) (< + l/t){-Q + cothsi coths2) 



D 



47-1 



D 



47-1 



1 



(3(3-^2 _ i/f2^ +cothsiCoths2 (3(-3- _ l/i2) + cothsi coths2 
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The following notation are used: 
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D 



(34) 



Q — (1 — 7) [(1 + 7) cothsi coths2 — (1 — 7) — i^Q;A(cothsi + coths2)] , 
L»_ = tooth si + (l/t)coths2 + \t-l/t\, D+ = tcothsi + (l/t)coths2 + {t + l/t). 



(35) 



The part of the vertex that is independent of Za requires the ultraviolet regularization and hence is dependent on 
the regularization parameter. In fact we have put the regularization parameter to be equal to 1/ro. The corresponding 
Z-independent term in Eq.(|3^) equals to ajlix. One should also take into account the term that provides the correct 
normalization of the total wave function, see e.g. Ref . |2^ , 



3 AT 



-a. J de J J Jdridr2dr3'ipp{r2)jf,gir2,r3\ie) 

— 00 

X 7"^ ('^3 , ri I ie)j''ipp (r 1 )D(r2 , ri\ie) 



(36) 



Straightforward calculation of this contribution with logarithmic accuracy gives —aLjl-K. Note that there are no 
.Z-dependent logarithmic terms in (^6|). As a result the "normalization" contribution cancels out exactly the Z- 
independent term in Ay . In essence this cancellation is a direct consequence of the Ward identity. The leading term 
of Za-dependent part of Ay at Za <C 1 reads 



Ay = 



(37) 



Numerical calculation of integrals in Eq. (34) together with account of the normalization contribution ( p6[ ) gives the 
function Ay exactly in Za. This function is plotted in Fig. ^ 
Equations (p^ and (|37|) give the final result for A at Za ^ 1 



A^A 



SE 



Ay = - 



a{Zaf 



15 _^ 



(38) 



The function A(Za) with account of all orders in Za is plotted in Fig. || 



V. LINEAR IN ZoL RADIATIVE CORRECTIONS 



In the leading a{Za) approximation, it is convenient to derive the non-logarithmic term of the radiative correction 
using the effective operator approach, see e. g. Ref. |28 . In this approach the corrections under discussion coincide 
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FIG. 4: The function Av calculated in all orders in Za. Value of the function is given in units 




FIG. 5: The function A calculated m all orders in Za. Value of the function is given in units 



with those for the cr • p structure in the amphtude of low-energy forward scattering. Diagrams for the scattering 
ampUtude are shown in Fig.^. Below we double the contribution of each diagram because of permutation of Z- 
boson and Coulomb lines. In these amplitudes the momentum |p| of the external electron is small compared to the 
momentum of the Coulomb quantum, k ^ m. Renormalization procedures for the self-energy operator and for the 
photon vertex operator are standard ones. For the Z-boson vertex we have to set the vertex correction equal to zero at 
k = because as a reference point we use result |13i th at is normalized at fc = 0. To minimize the infrared problems 
we perform calculations in the Fried- Yennie gauge |p3| where the photon propagator is of the form 



DM = 



(39) 



This gauge does not require an introduction of the photon mass for infrared regularization. In this gauge the renor- 
malized self-energy operator has the form 



A'Kni^{l — p) \ 1 — /O 



In p 



(40) 



where p = 1 — j? jw? . Calculating the diagram FigJ^(a) with this self-energy operator we obtain the following 
contribution to the function B defined in Eq. (13) 



Ba = -Q;(Za) . 



(41) 
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abed 



FIG. 6: Self-energy and vertex radtattve corrections to the forward scattering amplitude. The zigzag line denotes Z-boson and 
the dashed line denotes photon. The cross denotes the nucleus. Contribution of each diagram must be doubled because of 
permutation of Z-boson and Coulomb lines. 



A straightforward calculation leads to the following expression for the renormalized correction to the photon vertex 
operator 



A„ = £ |-37„ ^ dy (^ln[l + Ty(l - y)] + + dy 



dx 



(42) 



dy 







"'^ (1 — x)xdx 







2b(pi + m)ja{P2 + m)b 



Here the incoming electron has the momentum p2 and outgoing electron the momentum pi = P2 + k, 

2p ■ h 

b ^ [pi{l - y) + P2y\ , T ^x + Ty{l~xy) + \y{l-x), r = — , A = —. 

Performing calculations with this vertex function we find the contribution of the diagram Fig.|^(b) 

The renormilized correction to the Z-boson vertex reads 



o^iZa) ( -l-^ln2 



L. = £ {-7. f dy (3 ln[l + ryU - y)] 4) + ^ / 



^ dx 



[b'^la + 267qS) + 2x (jaP2b + bpi-ia - haPl " P27ab + m{irfa - lab) 



dy 



""^ (1 — x)xdx 



25(pi + m)^a{P2 ~ m)b 



+ (-2 + x)^^{pi + m)^a{P2 - "^)7M] + 
Calculation with this vertex function gives the following contribution of the diagram Fig.^(c) 

Be = a{Za) 



75- 



1 10, „ 

n 2 

4 3 



The last diagram, Fig.^(d), does not require any renormalization. A straightforward calculation gives 



Bd = a{Za) 1^- - + 4 In 2 
Altogether, the result for the non-logarithmic term in Eq.(0) reads 



B^Ba+Bb + Be + Bd^ -a{Za) 
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2 In 2 



Thus, according to Eqs. (Q), (p8|), and ( ^ ) total relative correction in the leading Za-approximation reads 



5ef = 



{Za) ( Y^ + 21n2 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 



The curve corresponding to this formula is shown in Fig.|^ by the dotted line pncl. By the dashed line pnc2 we show 
the correction 5f,f calculated with the exact in Za function A. The solid line pnc in the same Fig.|^ shows the total 
radiative correction to the PNC effect that includes both S^f and 5b (Eq.(||), Fig.|l|(b)). The leading unaccounted 
contribution in 5ef is of the order of ~ Z'^a^ /tt. For Cs (Z=55) this gives about 5% uncertainty in 5ef. 
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VI. RADIATIVE CORRECTIONS TO THE FINITE-NUCLEAR-SIZE EFFECT 



The radiative shift of the atomic energy levels (Lamb shift) depends on the finite nuclear size. This correction has 
a very similar structure to that of the PNC radiative correction since the effective sizes of the perturbation sources 
in both cases are much smaller than A^. The self-energy and the vertex corrections to the finite- nuclear-size effect 
(SEVENS) for si/2-state have been calculated earlier analytically in order a{Za) in Refs. |P9| , M . The corrections 
for lsi/2-, 2si/2-, and 2p]^/2"States have been calculated numerically exactly in Za in Refs . |p9[ |33|, . However, the 
structure of higher in Za corrections and their logarithmic dependence on the nuclear size has not been understood. 
We have applied our approach to the SEVENS problem and found the following expression for the Si/2-state relative 
correction in the leading Za-approximation 



A,, = -c 



(Za)(^-41n2 



(Zar 



^ - ^ ) ln(6Ac/ro) 



(49) 



Linear in Za term agrees with results of Refs. The correction given by Eq. ( ^ ) is shown in Eig. ^ by the 

dotted line fsl. By the dashed line fs2 we show the correction A^ calculated with the exact in Za function A. Results 
of the computations [|l^, |2^ for Is and 2s states are shown by circles and diamonds, respectively. The agreement is 
excellent. 

Note that logarithmic terms (ln(6A(7/ro)) in Sef and A^. coincide. This statement is valid in all orders in Za. 
Moreover, the logarithmic term in the SEVENS correction Ap for pi/2 state is also equal to that in 5ef and A^. The 
reason for this equality is very simple. The logarithmic terms come from small distances (r << Ac) where the electron 
mass can be neglected. When the mass is neglected the relative matrix elements for the PNC radiative correction and 
for SEVENS are equal. 




FIG. 7: Relative radiative corrections (%) for the PNC and for the finite-nuclear-size effects versus the nuclear charge Z. The 
dotted line pncl shows the correction Sef (Fig^(e,f)) when both logarithmic and nonlogarithmic terms are calculated in leading 
in Za orders, see The same correction, but with logarithmic term calculated exactly m Za, is shown by the dashed line 

pnc2. The solid line pnc shows the total radiative correction to the PNC effect that includes both Sef and Sb (Eq.^, Fig]^(b)). 
The dotted line fsl shows the finite nuclear size correction As when both logarithmic and nonlogarithmic terms are calculated 
in leading in Za orders, see Eq. (^^. The same correction, but with logarithmic term calculated exactly in Za, is shown by the 



loqarittir 

dashed line fs2. Results of computations of Aa for Is and 2s states JT^ , [g^/ are shown by circles and diamonds, respectively 



The correction Ap is calculated in our work . The result reads 



Ap 



1 \ , X (ZaY /15 



(50) 



Structure of this correction is qualitatively different from that of and A^, the expansion starts from a\n{l/Za) 
term, while in Sef and A^ it starts from Za^. This difference is due to the different infrared dynamics, see discussion 
in 11. 

It has been suggested in Ref. Q that the relation Sef = (A^ -t- Ap)/2 is vahd. Our results (|8|), (|^), and (|5^) 
clearly disagree with this relation. Although, due to accidental compensations, the relation is more or less valid 
numerically around Z « 57, we strongly insist that there is no valid justification for the relation. The "derivation" in 
Ref. [ ^2[ is based on an assumption that there is a gauge in which the vertex contributions to Sef, A^, and Ap vanish 
simultaneously. This assumption is obviously wrong, and this is why the above relation is wrong. 
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VII. CONCLUSION 

Now we can perform a consistent analysis of the experimental data on atomic parity violation since all the contribu- 
tions are known. In our analysis for Cs we include the theoretical value of the PNC amplitude from Refs. 0, |[ |[ |Tc| ] 

EpNC = 0.908(1 ± 0.005) 10"" ieaB{-QwlN), (51) 

as well as the —0.61% correction due to the Breit interaction jO| , the —0.85% radiative correction calculated in the 
present work, the +0.42% vacuum polarization correction [|l6[ |17[|, the —0.2% neutron skin correction the —0.08% 
correction due to the renormalization of Qw from the atomic momentum transfer q ^ 30MeV down to g = 
and the +0.04% contribution from the electron-electron weak interaction [T^ . The theoretical uncertainty from the 
Z'^c? j-K term unaccounted in the present calculation is about 0.05 — 0.1%. In Ref. |^ the ratio -Epatc//? has been 
measured with 0.35% accuracy. Here /3 is the vector transition polarizability. An analysis of the recent data on /3 has 
been performed in Ref.[|l0|. We use the value /3 = 26.99(5)a'^ obtained in [|l0|. Combining all these results we obtain 
the following value of the nuclear weak charge Qw at zero momentum transfer 

Cs: Qw = -72.81 ± (0.28)^. ± {QM)theor- (52) 

This value agrees with prediction of the standard model, Qw = —73.09 + 0.03, see Ref. We have used the neutron 
skin correction in our analysis. However, in our opinion, status of this correction is not quite clear because data on 
the neutron distribution used in Ref. [p5[ are not quite consistent with data on the neutron distributions obtained 
from proton scattering, see e. g. Ref. p6| | 

In the analysis for Tl we have included the theoretical value of the PNC amplitude from Refs. as well as the 
—0.88% correction due to the Breit interaction the —1.48% radiative correction calculated in the present work, 
the +0.90% vacuum polarization correction jl^, the —0.2% neutron skin correction, the —0.08% correction due to 
the renormalization of Qw from the atomic momentum transfer q ^ 30MeV down to q = |Q, and the +0.01% 
contribution from the electron-electron weak interaction |l7| . Using these theoretical results we obtain from the data 
Q the following value of the nuclear weak charge Qw at zero momentum transfer 

Tl: Qv^ = -116.8 ±(1.2)e,±(3.4)tfteor. (53) 

This agrees with prediction of the standard model, Qw = —116.7+ 0.1, see Ref. Q. Values of the weak charge Qw 
obtained in the present work differ a little from that reported in [|8|. There are two reasons for the difference, a) in 
the present work we have used a slightly different value for Cs vector transition polarizability, this value is probably 
more accurate, b) in the present work calculation of the logarithmic term in the radiative correction is performed 
exactly in Za, while in jlj] we have used only the leading term of Za-expansion. 

Concluding, we have calculated the radiative corrections to the effect of atomic parity violation, the corrections 
are enhanced by the collective electric field of the nucleus . This calculation has allowed us to perform a consistent 
analysis of the experimental data on the nuclear weak charge. Agreement with the standard model is within 0.6(t. 

We would like to thank V. A. Dzuba for helpful discussions. A.I.M gratefully acknowledge the School of Physics at 
the University of New South Wales for warm hospitality and financial support during a visit. 
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